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The phase diagram of a water solution of DNA and oppositely charged spherical macroions is 
studied. DNA winds around spheres to form beads-on-a-string complexes resembling the chromatin 
10 nm fiber. At small enough concentration of spheres these "artificial chromatin" complexes are 
negative, while at large enough concentrations of spheres the charge of DNA is inverted by the 
adsorbed spheres. Charges of complexes stabilize their solutions. In the plane of concentrations of 
DNA and spheres the phases with positive and negative complexes are separated by another phase, 
which contains the condensate of neutral DNA-spheres complexes. Thus when the concentration of 
spheres grows, DNA-spheres complexes experience condensation and resolubilization (or reentrant 
condensation). Phenomenological theory of the phase diagram of reentrant condensation and charge 
inversion is suggested. Parameters of this theory are calculated by microscopic theory. It is shown 
that an important part of the effect of a monovalent salt on the phase diagram can be described by 
the nontrivial renormalization of the effective linear charge density of DNA wound around a sphere, 
due to the Onsager-Manning condensation. We argue that our phenomenological phase diagram 
or reentrant condensation is generic to a large class of strongly asymmetric electrolytes. Possible 
implication of these results for the natural chromatin are discussed. 

I. INTRODUCTION and purification. 



In the chromatin a long negative DNA double he- 
lix winds around a positive histoiajs octamer to form a 
complex known as the nucleosomdil. Presumably due to 
their Coulomb repulsion nucleosomes position themselves 
equidistantly along the DNA helix, forming a periodic 
necklace or a beads-on-a-string structure, which is also 
called a 10 nm fiber. The fact that the self assembly 
of this structure is very sensitive to the salt concentra- 
tion shows that electrostatic forces are important. In 
proper conditions, the 10 nm fiber self assembles (con- 
denses) into the 30 nm chromatin fibeij- which is the ma- 
jor building material of a chromosomeu. It is interesting 
to understand whether electrostatic forces alone are suf- 
ficient to form the 10 nm fiber and determine the range 
of its stability. This is one of the motivations to study 
theoretically and experimentally a model of "artificial 
chromatin", where histone octamers (which themselves 
are in a complex equilibrium with a population of indi- 
vidual histones) are replaced by identical positive hard 
spheres (see Fig. |^). In experimental realizations of ^^a*- 
tificial Ghxomatin" , spheres can be colloidal particle^TQ, 
micellesBo or dendrimersQ. Other motivations are related 
to many applications of such complexes. One of them 
is the gene therapy. Inversion of charge of DNA by its 
complexation with positive spheres should facilitate its 
contact with a usually negative cell membrane making 
penetration into the cell more likely. Another applica- 
tion we want to mention here is the possibility of mak- 
ing DNA guided nanowireaa by complexation of metallic 
spherical colloids with DNA. When DNA is replaced by 
synthetic polyelectrolyte (PE), other industrial applica- 
tions emerge from pharmaceutics to protein extraction 
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FIG. 1. A beads-on-a-string structure of the complex of a 
long DNA helix with positive spheres. We assume that the 
number of turns, m, of DNA necessary to neutralize a sphere 
is large, m > 1. 

In the previous paper^l we analytically considered a 
worm-like negative PE (DNA) with the total charge —Q 
and uniformly charged positive spheres with charge q. 
(Below we use the word DNA instead of PE postponing 
generalizations to the end of the paper). We studied their 
complexation in a water solution containing a concentra- 
tion p of DNA molecules and a concentration s of spheres 
together with their monovalent counterions and a finite 
concentration of monovalent salt which provides Debye- 
Hiickel screening with the screening radius rg. In the 
Ref. ^ we dealt only with large enough concentrations 
s, neglecting entropy of free spheres. We have shown 
that, at large s, the state of the system is determined 
by the ratio s/si{p), where Si{p) — q/pQ is isoelectric 
concentration at which the total charge of DNA is equal 
to the total charge of spheres. At s < Si{p) DNA wind- 
ing around a sphere overscreens it so that the net charge 
of the sphere becomes negative. Negative spheres repel 
each other and form a periodic necklace. In the opposite 
case s > Si{p), more spheres bind to DNA than neces- 
sary to neutralize it, thus inverting the charge of DNA. 
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The origin of this counterintuitive phenomenon of charge 
inversion is the correlation induced attraction of a new 
sphere to a neutral complex: a new sphere approaching 
a neutral complex pushes away neighboring spheres, un- 
winds some of DNA from them and winds it on itself. In 
other words, a new sphere creates an oppositely charged 
image in the complex and binds to it. It is this additional 
attraction which leads to charge inversion. In this case, 
the net charge of each sphere is positive. Positive spheres 
repel each other and also form a periodic necklace. 

Far from the isoelectric point, s = Si{p), complexes 
have a beads-on-a-string structure (see Fig. |l|). They 
are strongly charged and repel each other in the solu- 
tion, so that solution is stable. At the isoelectric point, 
s = Si{p), the sphere net charge and the charge of the 
whole complex simultaneously cross zero. In the narrow 
vicinity of the isoelectric point, the charges of complexes 
(and their Coulomb repulsion) are so small that their 
short range attraction between complexes due to corre- 
lations of solenoids of DNA on different spheres (see Fig. 
H) is able to condense them. Complexes form large and 
weakly charged bundles (see Fig. ^ 
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FIG. 2. Cross section through the centers of two touch- 
ing spheres with worm-like (gray) PE wound around them. 
At the place where two spheres touch each other (the region 
bounded by the broken line) the density of PE doubles which 
in turn leads to a gain in the correlation energy of PE seg- 
ments wound around the spheres. 




FIG. 3. Almost neutral complexes condense into macro- 
scopic bundles. 

At a given concentration p of DNA and growing con- 
centration s of spheres, DNA molecules experience con- 
densation (aggregation) at some critical concentration 



s = Sc{p) below the isoelectric point s = Si(j)), and 
remain in condensed state until another concentration 
s = Sci{p) above the isoelectric point, where the inverted 
positive charge of the complex becomes so large that con- 
densate dissolves. The electrophoretic mobility changes 
sign at the isoelectric point s = Si{p) located between 
Sc{p) and Sd{p). Inversion of sign of electrophoretic mo- 
bility of complexes of DNA with dendrimers was recently 
reported in Ref. 0. 

In Ref. ^, we dealt only with large s and concentrated 
mainly on the structure and charge of free necklaces far 
from the condensation domain. In this paper, we calcu- 
late the dependence of the critical concentrations sdp) 
and Sd{p) for all {p, s) plane and study details of the con- 
densation domain. Our main result is shown by the phase 
diagram in Fig. |[ In this diagram, an aggregate of DNA- 
spheres complexes exists only in the region surrounded 
by the lines Sc{p) and Sd{p) (the two solid lines). The 
domain of large enough concentrations s, studied Ref. ^ 
is shown by gray. In the gray area condensation domain 
is narrow and encloses isoelectric line Si{p). We call this 
part of phase diagram the "neck" . At smaller s and p, 
the condensation region is found to be much wider so 
that we call it the "body" . The dash line corresponds 
to the values of s and p at which complexes are neutral. 
Notice that, in agreement with result of Ref. ||, this line 
is essentially the isoelectric line for the "neck", but it 
is high above the isoelectric line for the "body", where 
substantial fraction of spheres is free. 
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FIG. 4. Phase diagram of a water solution of long DNA 
molecules and positive spheres in the (p,s) plane. The dot- 
ted line corresponds to the isoelectric composition s = Si{p). 
The dashed line corresponds to the concentration of spheres 
s ~ so{p), where an isolated DNA-spheres complex is neutral. 
The solid lines Sc{p) and Sd{p) define the external boundary 
of the region of existence of macroscopic aggregates. Plus 
and minus are the signs of the charge of free DNA-sphere 
complexes and their aggregates above and below the dashed 
curve. The area of large s studied in previous papeitl is shown 
by gray. 

The same diagram is shown in linear scale on Fig. 0. 
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At large s and p, the phase diagram is extremely sim- 
ple and centered around the isoelectric line. When we 
zoom at smaller s and p, symmetry with respect to the 
isoelectric line disappears and we see the fine structure 
of the condensation domain. There are two light gray 
regions next to lines Sc{p) and Sd{p), where condensate 
coexists with free necklaces. Between them there is the 
dark gray region were practically all DNA and spheres 
are consumed by the condensate and concentration of 
free necklaces is exponentially small. We also study the 
change in the fraction of free necklaces and their charge 
when we cross all these regions, for example, by increas- 
ing s, while keeping p constant (see Fig. |^ below). One 
of interesting results is that the inversion of the necklace 
charge happens in a relatively narrow range of s. 




a) b) 
FIG. 5. a) Phase diagram of the mixture of long DNA he- 
lix and spheres in linear scale, b) Domain of small s and p 
is enlarged. On both plots, the dotted and dashed lines and 
plus and minus signs have the same meaning as in Fig. ^ 
The areas of coexistence of aggregates and free complexes are 
shown by light gray. In the dark gray region, all DNA helices 
condense. It is shown only in plot b) because it should be 
very narrow in plot a). One can see that the transitions to 
and from the complete condensation are abrupt at small p 
but become less steep with growing p. 

In Figs. ^ and || we did not try to show that the lines 
Sc{p) and Sd{p) should merge at extremely small p, where 
translational entropy of DNA molecules becomes impor- 
tant and as a result DNA complexes dissolve. These lines 
may also merge again when s and p are so large that the 
fraction of water in the mixture is significantly reduced. 
Then our phase diagram can also be redrawn as a close 
loop ternary miscibility diagram. 

Bundles of complexes formed in the gray area of phase 
diagram of Fig. 5 are almost neutral in the sense that 
their charge is proportional not to the number of aggre- 
gated complexes or volume of the bundle, but to the lin- 
ear size or surface area of bundle. The sign of the charge 
is however same as for free complexes. It fiips on the 
neutrality line where charge inversion takes place for free 
complexes. 

To derive these phase diagrams, we use two phe- 
nomenological parameters. The first one, e, is the bind- 
ing energy per sphere of the DNA-spheres complex (neck- 



lace) in the aggregate. The second one, sq, is the con- 
centration of spheres, which is in equilibrium with neu- 
tral DNA-spheres complexes. This concentration sep- 
arates domains of positive and negative complexes at 
p = 0. These two phenomenological parameters can be 
extracted from the experimental phase diagram. They 
can also be calculated from microscopic theory. In this 
paper we present comprehensive microscopic theory of 
these parameters and study the influence of salt on their 
value. We show below that, generally speaking, screen- 
ing by salt increases the width of the neck and, at strong 
screening, increases sq so that the "body" grows at the 
expense of the "neck" . 

It is known that the Qjisager-Manning condensation of 
monovalent counterionsQ renormalizes linear charge den- 
sity of free DNA helix from the bare value — ?/ to critical 
Onsager-Manning density —77c = — D/csT/e. Thus one 
can ask which of the two one should use to calculate 
charge of DNA molecule, Q, and the isoelectric point 
Si. We show that, because of the positive charge of the 
spheres surface, some DNA counterions are released into 
solution. As a result, the absolute value of the effective 
net linear density, rj* , of DNA (which almost neutralizes 
a sphere) decreases with decreasing screening length 
from the bare density rj to the Onsager-Manning critical 
density r]c according to the formula: 

Inf 7* / cl] 

- )' y {aeMvhc) > rs > Ao/2tt), (1) 

hi[AQ / Zira) 

where Aq = AnrjcR^/Q is the distance between DNA 
turns on the surface of a sphere if 77* — rjc and a is the 
radius of a DNA helix (see Fig. p. One can show that 
rj* = r]c if rs < Ao/2'k. 

It is the net charge density rj* of DNA which de- 
fines a renormalized isoelectric concentration s*{p) — 
{rj* /ri)si(p). Thus, when decreases, the phase diagram 
moves down together with s*(p) as shown in Fig. ^. As 
we mentioned above, screening also leads to some change 
of the width of condensation domain, but this effect is 
negligible compared to the shift of the isoelectric line. 
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FIG. 6. The shift of isoelectric line and the "neck" of phase 
diagram (gray) to lower s as a result of the addition of mono- 
valent salt. The points 1 and 2 are used in Sec. IV. 
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The phase diagram presented above and its phe- 
nomenological theory actually have a much broader ap- 
plications beyond the complexation of DNA with large 
spheres. It was actually discovered in experiments of 
complexation ofpOppositely charged proteins more than 
half century agoEj but have not get a theoretical expla- 
nation. Other applications are discussed in Sec. IV. 
Here, we would like to mention only the well known phe- 
nomenon of reentrant condensation of double-helix DNA 
in solutions with small size multivalent counterions with 
charge Z > 3 (Z-ions), such as spermine, a short positive 
PE with Z = A. At some critical concentration of^Z-ions, 
Sc, DNA abruptly condenses into large bundlesM. (We 
use for Z-ions the same notation as for spheres to em- 
phasize complete analogy at the phenomenological level) . 
Recently, it was discovered that at another much larger 
critical Z-ion concentratioru--|S|iT, DNA bundles dissolve 
again into free DNA heliceal3~ll3. For spermine at very 
small DNA concentrations, Sc = sdp ^ 0) = 0.025 mM 
and Sd = S(i(p ^ 0) = 150 mM, if the monovalent salt 
concentration is not large. 

A theory of reentrant condensation of DNA by sper- 
mine was given m Ref. 0. It is also based upon two phe- 
nomenological parameters e and sq , which have meaning 
of the binding energy of DNA per Z-ion in the aggregate 
and the concentration of Z-ions which is in equilibrium 
with ncutralized-by-Z-ions DNA. By comparison of Sc 
and Sd with theory, they were found to be £ ~ O.SfcsT 
and So — 3.2 mM. Ref. 17, however, dealt only with 
very small p. The phenomenological theory of this pa- 
per extends the theory of the spermine induced reentrant 
condensation to the whole plane (s,p). 

Although DNA-spheres complexes and DNA with sper- 
mine have the same phase diagram, there are important 
quantitative differences between them. This is because, 
as we will see in the next section, the concentration sq 
decreases exponentially with the binding energy of a Z- 
ion or a sphere with DNA. For large strongly charged 
spheres, the binding energy of a sphere with wrapping 
DNA can be so large that sq is extremely small. There- 
fore, for the case of "artificial chromatin" , one most likely 
deals with the "neck" of the phase diagram of the Fig. 
^ or, in linear scale, with Fig. Bk. This agrees with the 
phase diagram obtained in Ref. for DNA-large colloids 
complexation. On the other hand, for small Z-ions, the 
energy of their binding to DNA is much smaller and sq 
is an easily observable concentration. In this case, one 
should see Fig. |b or the left side ("the body") of the 
phase diagram of Fig. ^ instead. Indeed—the experimen- 
tal diagram obtained for the sperminelij looks like the 
left half of Fig. |. 

The paper is organized as follows. In Sec. II, we 
present a phenomenological theory of the complexation 
of DNA with spheres. Analytical formulae for the critical 
concentrations Sc{p) and Sd{p) are derived and details of 
the shape of the condensation and decondensation tran- 
sitions near sdp) and Sd{p) is discussed. In Sec. Ill 
we present microscopic theory for parameters e and sq. 



for artificial chromatin and study the role of screening 
by monovalent salt. Onsager-Manning condensation on 
DNA wrapping around a sphere is considered in Sec. IV. 
Only in Sec. IV, we return to the discussion of other ex- 
amples of application of the same phase diagram. In Sec. 
V, we discuss limits of applicability of our theory of reen- 
trant condensation. In the conclusion, we summarize our 
results and discuss the possible implication of this theory 
to the natural chromatin. 



II. PHENOMENOLOGICAL THEORY AND 
PHASE DIAGRAM OF REENTRANT 
CONDENSATION 



To begin with, let us first study a complex of a single 
DNA with N spheres. Its free energy can be written as 



f{N) = Q* '/2C + NE{N) 



(2) 



where Q* = qN - Q = Q{N/Ni - 1) is the net charge 
of the complex {q and Q are the charges of one sphere 
and a DNA molecule and Ni — Q/q is the number of 
spheres necessary to neutralize a DNA helix) . Except for 
a specific value of the sphere concentration s, the DNA- 
spheres complex is always charged and, therefore, has 
an extended shape to minimize its Coulomb self-energy. 
At a length scale larger than the average distance I be- 
tween the spheres along the complex, the complex can be 
considered as a charged cylinder with length L, radius I 
(see Fig. |l|) and linear charged density Q*/L. Then the 
capacitance C of the complex is: 



C 



DL 



21n(r,// + 1) 



(3) 



where D is the dielectric constant of water and monova- 
lent salt in solution is treated in the Debye-Hiickel ap- 
proximation with screening length r^. 

The first term in Eq. (|^) is the standard self-energy of 
a complex with net charge Q* and capacitance C. The 
second term is the correlation energy, which accounts for 
the discreteness of the spheres charge at the length scale 
smaller than I (in the beads-on-a-string structure, it is 
essentially the interaction of a sphere with the DNA coil 
wound around it). The correlation energy per sphere 
E{N) is negative and, in general, it is to be calculated 
from microscopic theory. In this section, it is assumed to 
be known and is used as an input parameter of the the- 
ory. Below, we show that the parameter sq mentioned in 
the introduction is directly related to E{N). 

Given a DNA concentration p, we want to calculate the 
value of the sphere concentration sdp) and Sd{p) where 
condensation and decondensation of DNA-spheres com- 
plexes happen. To do so, let us consider the system in a 
transition state where aggregates coexist with free com- 
plexes. Let X be the fraction of DNA in the aggregates, 
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the concentration of free DNA-spheres complexes in so- 
lution is then (1 — x)p. Correspondingly, the concentra- 
tion of free spheres in solution is s — xpNi — (1 — x)pN, 
where xpNi is the concentration of spheres consumed by 
the neutral macroscopic aggregates, and (1 — x)pN is the 
concentration of spheres bound to the free DNA-spheres 
complexes in the solution. The free energy per unit vol- 
ume of the system can then be written as: 

F{N,x;s,p) = {l-x)pf{N)+xpN,[EiN,) + e] + 

ksT [s - xpN, - (1 - x)pN] 1^ li^^WV^^ii^^)?^. 

e 

(4) 

The first term in Eq. is the free energy density of the 
free complexes with f{N) given by Eq. (|^). The second 
term is the free energy density of the aggregates, iV^e is 
the energy gained per complex by forming the aggregates 
(compared to a free neutral isolated DNA-spheres com- 
plex in solution). This energy gain originates from the 
correlation-induced short range attraction between com- 
plexes mentioned in the introduction (e is negative, see 
Fig. 1^). Like E{N), e is to be calculated from micro- 
scopic correlation theory, but in this section, it is consid- 
ered as another input parameter of the theory. The third 
term in Eq. (Q) is the free energy density of the con- 
centration s — xpNi — (1 — x)pN of left over free spheres 
in solution. This concentration is assumed to be small, 
so that the solution of free spheres is ideal with the nor- 
malizing volume vq. The translational entropy of DNA 
is neglected in zeroth order approximation. This is valid 
if the DNA helix is long enough and the concentration p 
of DNA is not too small. 

For given average concentrations p of DNA and s of 
spheres, the state of the system can be found by mini- 
mizing the free energy (||) with respect to the aggregate 
fraction x and the number of spheres N bound to a free 
DNA-spheres complex. This gives: 



fcBrin[(s - xpN, - (1 - x)pN)v^] = n,{N) + ^ (5) 

iv.. jqw^.-"'% w-iv)x 

{kBT\n[{s- xpN,-{l-x)pN)vo\- ^Jic{N,)) . (6) 

Eqs. and (|^) have very simple physical mean- 

ings. Eq. (||) equates the chemical potential of free 
spheres in solution (the left hand side) and the electro- 
chemical potential of spheres bound to the DNA-spheres 
complex (the right hand side). Here Q* /C is the aver- 
age electrostatic potential at the surface of the complex 
and Aic(iV) = d[NE{N)]/dN < is the contribution to 
the cherpjcal potential due to sphere correlations in the 
complexEl. 

Eq. (^ is the equilibrium condition for neutral com- 
plexes. The left hand side of Eq. (|^) is the binding 
energy per complex in the aggregates and the right hand 



side is the chemical potential of a free DNA-spheres com- 
plex in solution. The latter chemical potential is the sum 
of the complexes' self-energy (neglecting its translational 
entropy) plus the entropy of {Ni — N) spheres released 
into solution. (In the dissolved state, a neutral complex 
releases (Ni — N) spheres charging itself to the charge 
Q{N/N,-1).) 

Before progressing further, it should be mentioned 
that, both lengths L and I are, generally speaking, func- 
tions of the number of spheres in the complex, TV. A 
detail calculation on the dependence of L and I on N 
is given in Ref. ^, where the structure of complexes far 
from isoelectric point (deep in the plus and minus re- 
gion of the phase diagram of Fig. ^) was considered. In 
this paper, however, we concentrate on the phase dia- 
gram of the reentrant condensation. In the vicinity of 
the condensation region, the complexes are almost neu- 
tral. Therefore, in the above minimization of the free 
energy, L and / (and hence the capacitance C) are con- 
sidered as independent of N and equal to their value for 
a neutral complex. Indeed, if we consider explicitly the 
dependence of L, I and C on N, on the right hand side 
of Eq. (|^) we get an additional term of the order Q* /Q 
compared to the second term. It is shown below that 
Q* /Q is proportional to the square root of the ratio 



\eC/qQ\ < 1 



(7) 



where s is the condensate binding energy per sphere and 
qQ /C is the interaction energy of a sphere with the whole 
DNA molecule. Thus, the dependence of L, I and C on 
N can be ignored. The inequality (^ is justified because, 
as shown in Fig. |^, e involves only interactions with a 
small part of DNA molecule winding around a spheres. 
A more quantitative justification of this assumption is 
given in the next section where microscopic theory of 
DNA-spheres interactions in the complexes and the con- 
densate is given. 

If we define a sphere concentration sq as 



So = exp{~\fj.c{Ni)\/kBT)/vo 



Eq. (^ can be rewritten as 



xpNi - (1 - x)pN 
Sq 



N 



- 1 



qQ 
C 



(8) 



(9) 



From this equation, one can see that, when the concen- 
tration of free spheres, s — xpNi — (1 — x)pN, is greater 
than the concentration sq, N/Ni — 1 is positive indicat- 
ing a charge inversion effect. Because Coulomb interac- 
tion between spheres is much larger than thermal energy, 
|/ic(TV)| ^ ksT, So is an exponentially small concen- 
tration. The range of s where DNA-spheres complexes 
are overcharged, therefore, is easily accessible experimen- 
tally. One also sees that, without correlations, /ic ~ 0, 
So = I/wq and charge inversion can never be observed. 
This confirms that correlation is the driving force for 
charge inversion. 
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Using Eq. ^ and the assumption that the charge of a 
complex near condensation region is almost zero so that 
IJic{.N) ~ ^ciNi), we can rewrite Eq. (0) as 



e = 



N 



- 1 



qQ 
2C 



(10) 



Equations (||) and (|T^) make up the core of our the- 
ory. They contain only two phenomenological param- 
eters £ and So. The former describes the strength of 
short range attraction between complexes and the later 
describes the strength of sphere-DNA binding, both de- 
pend on the charge and size of spheres. Knowing them, 
one can easily calculate the fraction x of DNA in the 
condensate and the net charge Q* — Q{N/Ni — 1) of 
free complexes as functions of s and p. Thus, the whole 
phase diagram of the reentrant condensation can be con- 
structed. 

Let us consider the limit of small DNA concentration 
p. In this case, neglecting p inside the logarithmic func- 
tions, we can rewrite Eqs. (^) and ( |l0| ) as 



|£| kBTC^^^^ 



knT 



2qQ 



So 



(11) 

(12) 



Obviously, there are two solutions Sc and Sd for Eq. (12) 



ScM So exp T 



knT 



C 



(13) 



The solution with the minus sign, Sc < So, corresponds 
to the concentration of spheres at which DNA complexes 
start to condense forming aggregates. The other solution 
with the plus sign, Sd > so, corresponds to the concen- 
tration of spheres at which the aggregates dissolve again 
into free DNA-spheres complexes in solution. In Fig. ^ 
and Fig. the concentrations so, Sc and Sd are shown 
on the vertical axis. 

The meaning of the concentration sq is also transparent 
in this p — » limit. It corresponds to the sphere concen- 
tration at which a free complex is neutral [N — Ni). 

Let us now study the dependence of the threshold con- 
densation and decondensation concentrations Sc(p) and 
Sd(p) on the DNA concentration p. Because the bind- 
ing energy e is a constant, Eq. ( [lO|) shows that, in the 
regime of coexistence of aggregates and free complexes, 
TV is a constant with respect to p and s. Let us denote 
this constant by Nc for the condensation transition and 
by Nd for the decondensation transition: 



Nc4 ^nAIt ^mC/qQ 



(14) 



The deviations of Nc and Nd from Ni is related to the ef- 
fective charges of the DNA-spheres complexes when they 
are in equilibrium with their aggregates: 



Qld = Q{Nc.d/N, - 1) = TQVmc/qQ ■ (15) 

Together with the assumed inequality (^, Eq. ( p^ ) 
confirms that the net charge of a complex near conden- 
sation region is very small compare to the bare charge of 
the DNA. 

Eq. (||), in turn, leads to the conclusion that in 
the coexistence regime, the concentration of free spheres 
[s — xpNi — (1 — x)pNc,d] in solution is a constant inde- 
pendent of p. These constants, of course, are equal to 
the concentrations Sc and Sd obtained above for the limit 
p — > 0. Thus, for the condensation transition, 



s{p; x) - xpNi - (1 - x)pNc 



(16) 



and correspondingly, for the decondensation transition 



s(p; x) - xpN, - (1 - x)pNd = Sd ■ 



(17) 



Now one can easily find the threshold concentrations 
Sc{p) (where DNA-spheres complexes start to condense 
into bundles) and Sd{p) (where all the bundles dissolve 
again). Putting a; = in Eqs. ( p^ ) and ( [T^ ) one gets: 



Sc{p) = Sc+pNc 

Sd{p) = Sd +pNd 



(18) 
(19) 



Thus we have simple linear expressions for the condensa- 
tion and decondensation threshold sphere concentrations 
as functions of the DNA concentration. These functions 
are plotted in Fig. || by the solid lines. The same phase 
diagram is plotted in log-log scale in Fig. ^ where the 
lines distinguishing the coexistence and the complete con- 
densation regions are omitted for clarity. 

Putting a; = 1 in Eqs. ( |l6| ) and (^, one can easily cal- 
culate the two sphere concentrations between which the 
aggregates consume all the DNA-spheres complexes (we 
are talking about very long DNA neglecting its transla- 
tional entropy): 



s'M 

s'diP) 



Sc{p;x - 
Sd{p;x 



I) 



Sd 



-PN^, 

-pN,. 



(20) 
(21) 



One concludes that the width of the range of concentra- 
tion s where 100% DNA is aggregated, s'^{p) — s'^p) — 
Sd ~ Sc, is a constant. Remarkably, this region actually 
does not enclose the isoelectric line = pNi. In Fig. ||b, 
the complete condensation region is shaded dark gray. 

The sphere concentration so (p) at which a free complex 
is neutral can also be calculated in this model. Setting 
N = N, in Eq. (|) gives 



so{p) = So +pNi 



(22) 



As expected, sq{p) lies in between s' (p) and s'd{p) (com- 
pare Eq. (ll) with Eqs. (|o|) and (|lj)). In the phase dia- 
grams of Fig. ^ and ||, the concentration so{p) is plotted 
by the dash line. At large enough p, this concentration 
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is close to the isoelectric point. At small p, however, it 
saturates at the finite concentration sq. 

It is useful to study the relative width of the conden- 
sation region: 



As(p) sd{p) - Sc{p) srf - Sc + p{Nd - Nc) 



so{p) 
At small p 

so{p) 



so{p) 



So + pNi 



Sd - Sc. 
So 



= 2 sinh 



kRT 



C 



(23) 



(24) 



For large spheres, Coulomb interactions are much larger 
than the thermal energy, the argument of sinh function 
is large, so that this relative width is exponentially large. 
It is shown as the "body" of the phase diagram (Fig. |j). 
On the other hand at large p 



Asjp) _ Nd~Nc 
sq{p) 



N, 




< 1 



(25) 



The narrowing of relative width As{p)/so{p) with 
growing p from the "body" to the "neck" of the phase 
diagram is clearly seen in the log-log scale of the phase 
diagram (Fig. ^) . (Note that the "neck" encloses the iso- 
electric line as Eqs. ( |l8| ) and ( [l^ ) suggest.) The width of 
the condensation region, of course, always increases with 
p as shown in Fig. 0. 

We would like to note that, because of the exponential 
dependence of the concentration sq on the correlation 
chemical potential of spheres fic (see Eq. (^)), the ac- 
cessible range of the diagram is very sensitive to the size 
and charge of spheres. When spheres are large and highly 
charged, l/id/ZcsT is very large, sq is unrealistically small 
and one can access only the right side of this phase di- 
agram (the "neck"). For spheres of smaller size sq may 
be a more reasonable concentration and one can see also 
the left half of the phase diagram ( "the body" ) . 

Given the sphere and DNA concentrations, s and p, 
one can also calculate the fraction (1 — x) of DNA dis- 
solved in the solution. 



1 



1 - x{s) ~ < 







s - Sd - St 
Sr{Nd/N, - 1) 



1 



S < Scip), 

s'diP) < « < Sd{p), 

s'cip) < s < s'^ip), (26) 

s'diP) < ^ < ^d{p), 
S > Sd{p)- 



This result is plotted in Fig. ^ by the solid line. The 
slopes of x{s) in the condensation coexistence region and 
in the decondensation coexistence region are equal in 
magnitude. They are inversely proportional to the DNA 



concentration p. Therefore, at small p, the transitions 
from a; = to a; = 1 and vice versa are very abrupt. 
At large p, the coexistence region is wider. This fact is 
shown in the phase diagram of Fig. ||. The widths of the 
light gray regions corresponding to coexistence are zero 
at p and increase linearly with p. In the limit of 
large p, the function 1 — x(s) acquires a V-shape form 
between Sc{p) and Sdip)- 




FIG. 7. The fraction (1 — x) of free DNA in solution and 
the charge Q* (in units of Q) of a free complex as a function 
of the concentration s of spheres. 

Let us now consider how the charge Q* of a free com- 
plexes varies from negative at s < So{p) to positive at 
s > soip). When s < sdp), Eq. (^ shows that Q* grows 
logarithmically with increasing s. However, in the first 
coexistent region, Sc{p) < s < s'^{p), the net charge of the 
complex Q* equal Q* < and stays constant. In the sec- 
ond coexistence region, s'^{p) < s < Sd{p), the net charge 
is once more constant and but equal Q*d = —Ql > 0. 
When s > Sd{p), this net charge increases again logarith- 
mically with s. In Fig. 0, this net charge is plotted by 
the dashed hne. In the region s'^{p) < s < s'^{p), because 
the condensate consumes all the free complexes, we have 
to go beyond Eqs. (^ and ( p^ ) to find Q*. Namely, in 
this interval we should take into account translational en- 
tropy of free complexes. This gives a new equation which 
replaces Eq. (p^: 



^. - ksT-lnVl - x)pV„l (27) 



N.M ^ I I - 1 



where (1 — x)p is the concentration of free complexes in 
equilibrium with the condensate and Vq is the normaliz- 
ing volume for the complexes. 

Since we are interested in the variation of the net 
charge Q* = Q{N/Ni — 1) with respect to s, let us elim- 
inate X from the two Eqs. (||) and (p7|). This gives: 
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s — pNi = So exp 



qQ* _^N,vo ( N 



For large Ni (long DNA) , the second term in the above 
equation is exponentially small compared to the first 
term, thus one has: 



Q* 



ksTC 
q 



In 1 



So 



(29) 



Using Eqs. (|l^) and ( |l5|) one can verify that Q* matches 
Q* and at s = Sc{p) and s — Sd[p) respectively. 

It is interesting to note that the region where the com- 
plexes inverted charge is relatively narrow. Indeed, the 
slope at which Q* changes sign from negative to positive 
at so(p) is 



dQ* Q kBT2C 



ds 



So qQ 



(30) 



Thus the characteristic width Ss — 2Q*^/ [dQ* / ds) over 
which charge inversion happens is exponentially smaller 
than Sd — Sc- 



III. MICROSCOPIC THEORY 

In the previous sections, we presented a description 
of the correlation-driven reentrant condensation and ob- 
tained an universal phase diagram based on two phe- 
nomenological parameters e and sq. These two parame- 
ters depend on microscopic properties of the specific sys- 
tem considered. 

In this section we use a microscopic theory of the 
beads-on-string structure (Fig. |l|) for calculation of e. 
So and their dependence on monovalent salt concentra- 
tion. We consider large spheres and, therefore, effectively 
work in the "neck" of the phase diagram (large p case) . 

First, for simplicity, we start from the theory for poly- 
electrolyte with linear charge density 77 < rjc = ksT/e. 
The next section considers highly charge polyelectrolyte 
and the Onsager-Manning condensation which renormal- 
izes the net charge of DNA resulting in a shift of the 
isoelectric line. 

Because the DNA coil winding around a sphere almost 
neutralizes the sphere charge, the correlation chemical 
potential is essentially the self energy of a bare free 
sphere in solption which is almost totally eliminated in 
the complexu. Thus one has 



_ / -q^/2DR Ts > R, 
-q^rsl2DR^ Vs < R. 



(31) 



The parameter sq can be easily found with the help of 
Eqs. (H) and (|3l|). One sees that, at large Ts, the con- 
centration So is almost constant but when Ts decreases 
below i?. So increases exponentially with rg. 



To calculate the correlation attraction energy e, let us 
use Fig. When two spheres touch each other, the den- 
sity of the solenoid at the touching region doubles. This 
leads to a gain in the correlation energy of DNA turns 
(these correlations develop at distances of the order of A 
and should not be confused with the correlation between 
spheres which determines fic and develops at distance of 
the order of R) . 

The correlation energy per unit length of the DNA 
can be estimated as the interaction energy of the DNA 
segment with its stripe of background (sphere) positive 
charge. Thus it is —rj'^ ln{R/A)/ D for > R and 
—ri'^ln{rs/A)/D for R > > A. Correspondingly, when 
the density doubles {A halves), the gain in the correla- 
tion energy per unit length is —rf/D for > A. Similar 
effect takes place at Ts < A. In this case, each DNA inter- 
acts with the stripe of the width of of the sphere sur- 
face with energy ~qr]rs/2DR^. When two spheres touch 
each other, it interacts with the stripe on the other sphere 
as well, doubling the correlation energy. Thus in this 
case, the correlation energy gain is —qrirs/2DR^. Sim- 
ple geometrical calculation shows that the total length of 
DNA in the touching region (the region surrounded by 
broken line in Fig. ^ is R. Therefore, 



-Rrf/D 
-qrirs/2DR 



< A. 



(32) 



Using Eqs. ( |3l| ) and (p2|), we are now in a position 
to discuss in detail the change of the phase diagram 
with varying monovalent salt concentration (varying r^). 
Let us consider the following three regimes: > R, 
R> rs > A and Vs < A 

In the first regime, > R, the parameters so and e 
remain constant. The capacitance C of a complex has 
the form: 



C = L»L/21n(r,./i?) , > R , 



(33) 



where in Eq. (^) the average distance I between spheres 
has been replaced by 2R, because, near the condensa- 
tion region, complexes are almost neutral and spheres 
almost touch each other. Therefore, Eq. ( p3| ) and ( |lO| ) 
show that, in the coexistence regions, the complex charge 
Q{Nc.d/Ni — 1) increases logarithmically with decreasing 
Tg. Eq. ( [25I ) then shows that the relative width of the 
condensation region increases slowly in this regime: 



As{p) 2R-q 



1 



sq{p) 



q V\n{rs/R) 



< 1 



(34) 



Here, Eq. (||) and the equality Q/L = q/R were used. 
Recall that we assume the number of turns, to, of DNA 
around a sphere is large. Therefore Rrj/q ~ to^^ <C 1. 
Decreasing Ts, however, leads to the contraction of the 
condensation region in the limit p ^ 0. Indeed, Eq. ( p^ ) 
shows that at constant e, the concentration Sc increases 
and the concentration Sd decreases. Thus, overall, the 
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condensation region decreases linearly with decreasing rg 
at small p and increases logarithmically at large p. 

In the second regime, R > rg > A, The concentra- 
tion So increases as exp(— r^/i?) while e remains constant. 
The phase diagram moves upward with sq. Since < i?, 
the capacitance C is 



C = DLR/2rs 



< R 



(35) 



Thus, the relative width of the condensation region in- 
creases faster than Eq. (BJ): 



As{p) 2Rr] 



(36) 



(Notice that Eq. (|6|) matches Eq. (||) at - R). 
At the same time, at small p, the condensation region 
shrinks faster Sc^d oc exp(=prs/i?). Thus in this regime, 
the change in shape the phase diagram continues the 
trend in the first regime but with a faster pace. 

In the third regime, rg < A, the energy e also starts to 
decrease linearly with r^. Eqs. ( p^ ) and ( p^ therefore 
lead to the saturation in the change in ScA and Nc^d- The 
relative width of the condensation region in the "neck" 
saturates at: 



so{p) 



l2Rr] 



(37) 



(Notice that, because A - R'^rj/q, Eqs. (|6|) and (H) 
match each other at rg ~ A.) On the other hand, the 
concentration sq continues to increase exponentially with 
decreasing r^. One can say that in this region, the main 
effect of increasing salt concentration is the growth of the 
"body" of the phase diagram Fig. ^ at the expense of the 
"neck" . 

The evolution of the width of the "neck" as function 
of rg is shown on Fig. 0. 



As 
So 

,-1/2 \ 



r.\-i/2 




-1/2 







A R 

FIG. 8. The relative width of the condensation region of 
DNA-spheres complexes as a function of the screening length, 
m — q/rjR S> 1 is the number of turns of DNA neutralizing a 
sphere. 

Concluding this section, we would like to justify quan- 
titatively the assumption of Eq. (^ , which leads to con- 
clusion that the charge of complexes is small near con- 
densation region. Using Eq. (|3^) for e and Eqs. (|3|) and 
Eq. (Ba) for the complex capacitance C one gets: 



\e\C 
qQ 



{Ri^/qYlHrg/R) Tg > R, 
{Rrj/qfR/rg R>r, 
(Rrj/q) A>r, 



>A^ R^Tj/q, 



(38) 



One sees that in all cases, the ratio \eC/qQ\ is at most 
Rrj/q = m^^ where m is the number of turns of DNA 
coil on a sphere which we assume to be large. Thus the 
inequality (|^) is justified. 



IV. ONSAGER-MANNING CONDENSATION 



The theory presented so far is applicable to the com- 
plexation of polyelectrolyte with linear charge density 
which does not exceed the so called Onsager-Manning 
critical charge density rjc = DkBT/e. DNA helix is, how- 
ever, a highly charged polyelectrolyte with the bare neg- 
ative charge density —j] ~ —A.2r]c. Condensation of pos- 
itive monovalent ions on DNA, therefore, must be taken 
into account. This section deals with the problem of 
counterion condensation. 

For a free DNA helix, the Onsager-Manning conden- 
sation results in the net charge —rjc- However, the net 
linear charge density of DNA winding around a charged 
sphere, —rj*, may differ from —rjc due to the repulsion of 
counterions from the positive charge of the sphere. This 
phenomenon is called counterion releasal3. We show be- 
low that, for DNA coil which neutralizes a sphere, rj* is 
given by Eq. (P, so that i] > i]* > ric- Similar to L and 
I, the net charge density of DNA, ?/*, can be considered 
as independent on s, p and A^, if one is interested only in 
the reentrant condensation domain of the phase diagram 
where the charge of complexes are small. Therefore, the 
theory used in previous section remains valid provided 
one replaces Q hy Q* = Qr]*/r], Ni by Q*/q and Si by 
pQ*/q. 

Remarkably, because the relative width of the conden- 
sation domain studied in Sec. Ill is small with or without 
screening, when Onsager-Manning condensation is taken 
into account the main effect of screening in the phase 
diagram is the decrease of Q and, correspondingly, the 
isoelectric concentration s^. The whole right part of dia- 
gram of Fig. m shifts in the direction of lower s together 
with the renormalized Si, as it is shown in Fig. To 
understand the implications of such shift, let us consider 
a solution I, which is represented by the point {s,p) be- 
low the curve Sc{p) at very large rg, when 77* is equal to 
the bare DNA charge density 77. If we add more mono- 
valent salt and reduce Vg, both curves Sc{p) and Sd{p) 
move down following the move of the isoelectric line. As 
the result, our point 1 can be found inside the domain 
where complexes aggregate. Thus condensation of nega- 
tive (underscreened) complexes of DNA and spheres can 
be induced by addition of salt. 
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On the other hand, if we start from the solution rep- 
resented by the point 2, which at large Ts is above the 
decondensation line Sd{p), the result is completely dif- 
ferent. In this case, addition of salt can not condense 
DNA-spheres complexes. The origin of this asymme- 
try is that we are considering large enough spheres with 
surface charge density smaller than that of DNA helix, 
so that only monovalent cations experience the Onsager- 
Manning condensation. 

It should be noted here that all the discussion above 
about the change in the isoelectric line is relevant for 
the "neck" of the phase diagram only. The expansion of 
the "body" toward higher p and s with decreasing at 
rs < R mentioned in the previous section is not affected 
by the Onsager-Manning condensation. 

Let us now derive Eq. (^) for the net linear charge den- 
sity of the DNA due to Onsager-Manning condensation. 
Let us denote by ci the concentration of monovalent salt 
in the solution and by the concentration of monova- 
lent ions condensed on the DNA surface. As we know, 
for a free DNA in solution, counterions condense on DNA 
making its net charge — r/c- The balance of the chemical 
potential of counterions in solution and those condensed 
on the DNA helices, then reads 



c, 2e?7r r, 
kB-L In — = — — In — , 
Cl D a 



(39) 



where the left hand side is the entropy lost and the right 
hand side is the energy gained when a counterion con- 
denses on the DNA surface, a is the radius of the DNA 
double helix (not to be confused with the radius I of 
the complex, which, for almost neutral complexes where 
spheres almost touch each other, is roughly the radius of 
one sphere, 2 ~ i? » a) . 

In the DNA-sphere complex, the positive charge of the 
sphere reduces the amount of the counterions condensed 
on the DNA making the net linear charge of DNA equal 
—77*, instead of —rjc- The equilibrium condition for the 
monovalent counterions now reads 



fc^Tln — = 

Cl 



2e77* , A/27r 
In 



D 



sV^O). 



(40) 



Here A/2tt — 2R^t]* /Q plays the role of the local screen- 
ing length at the DNA surface [A is the distance between 
subsequent turns of DNA around a sphere (see Fig. ^). 
If A/2tt ^ Ts then ^/27r should be replaced by r^.. The 
average electrostatic potential at the surface of a DNA- 
sphere complex is V'(O) = [2Q*{N/N, - l)/DL\ \n{rs/R). 
(Here and below, bearing in mind that we are working in 
the vicinity of isoelectric point we replaced Ihy R). 
Combining Eqs. (|9|) and (^), one gets 



(41) 



rir In — = 71* In ^— ( 1 | In — 

' a ' 2™ L \N, R 



Near the condensation domain, the charge of a complex 
{N/Ni — 1) is very small and the second term in the right 



hand side in Eq. (|4l]) is negligible. (Indeed, according to 
eq. d), Q*{N/N, - 1)/L ^ ksTD/q « keTD/e = r?„ 
and because R ^ a, the second term is much smaller 
than rjc ln(rs/a)). Using A/2t: = 2rj*R^ /Q and substitut- 
ing 77* by rjc in ln(A/27ra), one arrives at A^-independent 
net linear charge density for DNA given byEq. (|l|). This 
net charge decreases with decreasing (increasing salt 
concentration). When rg reaches A/2tt, the net charge 
density ry* reaches rjc- When is even smaller, A/2tt is 
replaced by and 77* saturates at the value of rjc- 

The above derivation of 77* is similar to that of Ref. 
pO| where the problem of adsorption of DNA (with its 
Onsager-Manning condensed counterions) on a charged 
surface was considered. However, the obtained net charge 
(Eq. (|l])) is different with that of Eq. (7) of Ref. This 
is because Ref. |2^ considers the case when DNA strongly 
overcharges the surface, while this paper concerns with 
almost neutral DNA-spheres complexes. 

In the final paragraph of this section, we would like 
to discuss the condensation of negative monovalent ions 
(colons of DNA) on the spheres. This happens if the 
surface charge density of the sphere is4arge enough and 
the screening radius is not very largeES. Near the con- 
densation domain, however, DNA coil almost neutralizes 
a sphere eliminating its Coulomb potential at distances 
larger than A/2tt. Therefore, most of colons are released 
from the sphere. There is still a small amount of them 
condensed in the middle between two DNA turns. How- 
ever, the total charge of these colons is much smaller than 
the total charge of DNA counterions condensed on the 
DNA because, for A':^ a, DNA has a much larger bare 
surface charge density. Thus these colons are not impor- 
tant in our calculation of the isoelectric point. This also 
means that our calculation of e is valid. On the other 
hand, /ic can be modified when sphere condensation is 
taken into account because the self energy of the sphere 
in this case is smaller than that given by Eq. (|3|) . How- 
ever, /ic affects only sq. For large spheres where Coulomb 
interaction is much larger than fc^T, we always are in the 
"neck" of the phase diagram and sq is irrelevant. There- 
fore, the effect of condensation of DNA colons on spheres 
is small. 



V. OTHER APPLICATION OF 
PHENOMENOLOGICAL THEORY OF 
REENTRANT CONDENSATION 



In Sec. II, we presented a theory of the phase diagram 
for complexation of a long DNA helix with oppositely 
charged spheres. It shows that correlation induced charge 
inversion leads to reentrant condensation. The phase di- 
agram is described by two parameters e and Sq- Be- 
side these two parameters, however, the phenomenologi- 
cal theory does not use any specific information about the 
structural and chemical properties of the DNA. There- 
fore, this theory is generic and can be used to describe 
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a broad range of systems experiencing reentrant conden- 
sation. In this section, we would like to discuss briefly a 
number of other systems of large opposite charges which 
can be described by the same phase diagram of Fig. ^. 

Spheres can be replaced by particles of any shape which 
we for brevity call Z-ions. There is no need for them to 
be rigid. For example, they can be star micelles or just 
oppositely charged shorter polyelectrolytes (PE), more 
flexible than DNA, so that DNA can be considered rigid. 
Still they will complex with DNA, for example, winding 
around it, repel each other on DNA and form a kind of 
necklace. Phenomenological parameters sq, e can still 
be introduced (but, of course, microscopic calculation of 
these parameters is different from that of Sec. Ill) and 
the phase diagram should look similar to Fig. ^ or Fig. 

i 

Let us discuss the limit when the Z-ion size is smaller 
than double helix diameter so that they form two di- 
mensional correlated liquid on the surface of DNA. A 
well studied example of this limit is the reentrant con- 
densation of DNA in the presence of spermine already 
mentioned in Introduction. Microscopically, it is differ- 
ent from the beads-on-string systems, because in this 
case DNA helices can be considered as rigid cylinders. 
Condensation of DNA, of course, requires attraction be- 
tween like-charged cylinders. This attraction is related 
to the fact that Z-ions, at the surface of DNA cylinders, 
strongly repel each other and form a two-dimensional 
strongly correlated liquid. At the place of the con- 
tact where the correlated liquids of two touching DNA 
merge, the two-dimensional concenJcmiion of Z-ions dou- 
bles. This leads to an energy gainEil EEl because the corre- 
lation energy of a strongly correlated liquid per Z-ion is 
known to be negative and increases in absolute value with 
increasing Z-ion concentration. In other words, two DNA 
helices with strongly correlated liquid of Z-ions on their 
surfaces experience a correlation-induced short range at- 
traction. [— in 

It was predictedE3E3 that in a very weak DNA solu- 
tion, the net charge of a DNA helix inverts sign at a 
critical concentration s — sq- It has the meaning of con- 
centration of Z-ions in solution which is in equilibrium 
with the strongly correlated liquid of Z-ions at the sur- 
face of DNA (concentration of "saturated vapour" above 
strongly correlated liquid ). Physics of this inversion of 
charge is also related to correlations: when a new Z-ion 
approaches already neutralized DNA it forms an image 
in strongly correlated liquid, which attracts it. 

Combining correlation induced charge inversion and 
short range attraction, Rcf. [ij offered an explanation for 
the origin of the reentrant transition at small concen- 
tration of DNA. At s < So charges of two helices are 
negative. Their absolute values decrease with increasing 
s until at s = Sc, where Coulomb repulsion looses to the 
correlation attraction and DNA condenses in bundles. 
At s > So the net charge of DNA becomes positive and 
grows with increasing s. At s = s^, the Coulomb repul- 
sion wins over correlation attraction and DNA bundles 



dissolve. For p ^ Eqs. (|lj), ( [T^ ) were derived already 
in Ref. What is done in phenomenological theory of 
Sec. II can be considered as extension of the theory of 
Ref. |l^ to arbitrary p. 

Applicability of the phase diagram of Fig. ^ of course, 
is not limited to DNA. For all the broad spectrum of Z- 
ions, long DNA helices, in turn, can be replaced by other 
long, strongly charged PE. Experiment on complexation 
in mixture of micelles and oppositely charged PEEI is a 
good example. In this experiment the total charge of 
micelles was varied by changing the concentration of the 
cationic lipid in the solution. At some critical point the 
electrophoretic mobility of complexes changes sign. In 
agreement with the above theory measurements of dy- 
namic light scattering and turbidity (coefficient of light 
scattering) show that complexes condense in bundles in 
the same vicinity of the point where mobility crosses over 
between two almost constant positive and negative val- 
ues. 

An interesting new case is that of a weakly charged 
(one charge per 1// ^ 1 monomers) flexible PE. Let 
us consider its complexation with small strongly charged 
spheres. In this case, when PE winds around a sphere, 
the entropy of PE chain prevents its collapse to the sur- 
face of the sphere. As a result PE builds around the 
sphere layers of blobs with changing size, similar to the 
one discussed in Ref. The effective radius of a small 
sphere can be renormalized to a larger value. The rest of 
theory of complexation is similar to Ref. |^ and this paper 
so that our phase diagram works in this case, as well. 

Complexation of a long PE with much shorter oppo- 
sitely charged PE in a water solution is important class of 
such problems. Let us assume that the long PE is flexible 
while the shorter one is rigid and stronger charged. Then 
the long one sequentially wraps around molecules of the 
short one so that rods-on-a-string system resembling Fig. 
is formed. Then the only change to be made in the the- 
ory of Ref. H and this paper is to replace self energy of 
a charged sphere by self energy of a charged rod. The 
case when the short polymer is flexible (but still stronger 
charged) on the first glance seems to be more compli- 
cated. However, away from isoelectric point short PE 
is either underscreened or overscreened by the wrapping 
long one, so that it has the rode-like shape. Similarly the 
whole complex is charged, stretched and always resem- 
bles Fig. |l|. Finally, if both PE are weekly charged but 
even the short one is long enough, so that its charge is 
large, they still form a complex resembling Fig. ^ with 
some renormalizations of geometrical parameters related 
to a hierarchy of blobs. 

A long DNA helices or long PE molecules can be also 
replaced by any large macroions, for example, a large 
colloidal particle, which adsorbs smaller Z-ions. With 
appropriate correction for the expression of the capaci- 
tance C of the complex which, in this case, has spherical 
instead of cylindrical shape, our theory is applicable as 
well. An interesting practical example of such system 
is solution of positive (latex) spheres with short DNA 
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oligomers with, say, Sjpr 16 bases, which are adsorbed at 
the surface of spheresdu. Modified DNA can be deUvered 
by such colloids as a drug. Therefore question of stability 
of such solutions at a given concentration of DNA is very 
important. 

It is interesting that in this case latex spheres play the 
role of long DNA helix in the theory of our paper, while 
short DNA oligomers play the role of spheres or Z-ions. 
Indeed, short DNA oligomers form two-dimensional cor- 
related liquid on the surface of colloid and can overscreen 
it. Simultaneously this liquid provides attraction of al- 
most neutralized colloids leading to reentrant condensa- 
tion around isoelectric point. Reentrant condensation 
and change of the sign of electrophoretic mobilitsi in these 
systems were carefully studied experimentallyS'cl and the 
results seem to agree with our theory. Electrophoretic 
mobility of free complexes changes sign in the very nar- 
row vicinity of isoelectric point where solution is unstable 
and easily coagulates. 

Condensation of colloids with multivalent Z-ioos. was 
recently observed in Monte-Carlo experimentsLj al- 
though the concentration of .Z-ions are not large enough 
to observe resolubilization. 

Another example of application of our phase diagram 
to almost spherical macroions is the solution of nucleo- 
somes with spermineli3. In the nucleosome the positive 
histone octamer (charge 160e) is strongly overcharged 
by DNA with charge — 292e. Adsorption of spermine 
cations on the surface of nucleosomes leads to nucleo- 
some attraction and aggregation. Phase diagram of this 

fstem looks like the body part of phase diagram of Fig. 
Electrophoretic mobility was measured for aggregates 
of nucleosomes. With increasing concentration of sper- 
mine it changes sign from negative to positive roughly 
in the middle between condensation and decondensation 
curves Sc{p) and Sc{p) (in logarithmic scale). 

VI. LIMITATIONS OF REENTRANT 
CONDENSATION. 

In the previous section, we saw that the combination 
of reentrant condensation and charge inversion is a very 
general phenomenon. Now we would like to understand 
whether this phenomenon takes place for any strongly 
charged and strongly interacting electrolyte. We start 
from a theoretical example of a system, which does not 
show this phenomenon. This will help us to formulate 
limitations for reentrant condensation and charge inver- 
sion. 

Let us consider complexation of an absolutely rigid 
cylinde&rlike negative PE with small and strongly charged 
Z-ioniE3. Let us assume that linear charge density of PE 
is 77, radius of the cylinder is a and it is larger than radius 
of Z-ion. We assume also that PE is not very strongly 
charged so that Z/ar] » 1. In this case, isolated com- 
plexes of a cylinder with adsorbed Z-ions resemble one- 
dimensional Wigner crystal. Two such one-dimensional 



Wigner crystals strongly attract each other. The deep 
minimum of Coulomb energy is always provided by the 
condensed state, where cylinders form densely packed 
background on the top of which Z-ions form a three- 
dimensional strongly correlated liquid with average dis- 
tance between Z-ions much larger than a (they actually 
sit in pores of the background). The binding energy of 
such strongly correlated liquid is much larger than for 
free one-dimensional complexes. Therefore, at large con- 
centrations of PE and Z-ions, this system is always in 
condensed state and there is no way to see reentrant con- 
densation. 



This (not very realistic) example emphasizes nontriv- 
ial nature of reentrant condensation and charge inversion 
in the systems of DNA with large spheres or with sper- 
mine (and other similar systems discussed above). The 
destruction of aggregates by Coulomb repulsion of com- 
plexes at small distance from neutrality line (reentrant 
condensation) is possible only if the condensate bind- 
ing energy per Z-ion, e, is smaller than the Coulomb 
interaction of a Z-ion with macroion, qQ/C. Inequal- 
ity e/{qQ/C) ^ 1 was used in Sec. II to simplify cal- 
culations. If it is violated, the net charge Q* ^ is not 
small, the neck widens and can occupy the whole plane 
at large s and p eliminating condensation and deconden- 
sation transitions. Of course, the requirement that inter- 
action between Z-ion and DNA is large (|/ic 3> fesTj), 
so that charge inversion can happen, is also necessary for 
an observation of reentrant condensation. 



Small values of e are realized for DNA with large 
spheres or with spermine because in both cases charges 
which form complexes very strongly screen each other. 
For DNA with large spheres charge of a sphere is well 
screened by coiling DNA so that spheres stick together 
only at the place they touch each other (Fig. ||) and the 
energy e is very small. In the case of DNA with sper- 
mine, Z-ions so strongly screen DNA that attraction is 
developed only in the small area where two DNA touch 
and the energy e is relatively small. 

Between systems of DNA with large spheres and DNA 
with spermine there are many others with intermediate 
size and charge of spheres. One can easily imagine an 
intermediate system, where on one side spheres are not 
large and strongly charged enough that DNA behaves as 
a flexible PE, and on the other side they are strongly 
enough charged so that they form something like one- 
dimensional Wigner crystal on its surface. In such a 
case, mutual screening of DNA and spheres is not as 
complete as in the case of DNA with large spheres or 
with spermine (see Fig. ^ or Ref. |3^). Therefore, in- 
equality \e\/{qQ/C) ^ 1 can be violated and reentrant 
condensation may not materialize. 
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FIG. 9. Complexation of DNA and spheres with interme- 
diate value of charge and size. Screening is not complete as 
in the case of large spheres or multivalent ions. 

Let us give an example of a different, well-studied sys- 
tem where this happens. Consider a solution cationic 
lipid membranes and long double helix DNA. Each mem- 
brane, in principle, may complex with many DNA rods 
which lie on it equidistantly parallel to each other, form- 
ing a kind of one-dimensional Wigner crystal. (In such a 
complex DNA plays the role of Z-ion) These complexes 
may be undercharged and overcharged. They can also 
condense in three-dimensional aggregates, where layers 
of equidistantly positioned parallel DNA rods alternate 
with membranesEll. Such aggregates have potential ap- 
plications in gene therapy as non-viral gene carriers. In 
spite of liquid nature of membranes, which permits some 
redistribution of lipid charge inside membrane to screen 
DNA, attraction between different complexes in this crys- 
tal is so strong that situation for a reentrant transition 
in this case is marginal {\e\C/qQ ~ 1). It seems from 
experimental data that indeed aggregates exist even far 
from the isoelectric pointEJ. 

For the colloidal spheres covered by short DNA men- 
tioned in Sec. V, binding is related to the small area of 
contact and energy e is small again. Let us, instead of 
short DNA, consider small positive colloidal particles. If 
their size and charge are much smaller than that of larger 
negative ones we still get reentrant condensation like with 
the short DNA. But what happens when their charges 
and sizes become comparable is not obvious. Therefore 
it is interesting to consider a simpler case of a solution 
of positive and negative spherical macroions with same 
absolute values of charges, g, same radius, R and with 
concentrations c+ and c_ respectively. Let us talk about 
large enough c+ and c_, when Coulomb energy of in- 
teraction of two touching spheres, E = q^/2DR, com- 
pletely dominates the entropy. It is obvious then that at 
c+ = c_ such system tends to self assemble into NaCl-like 
crystal and gain the Madelung energy —1.7AE per neg- 
ative sphere. Away from the isoelectric point situation 
is less trivial. Let us consider for simplicity a solution 
with c+ = 2c_ (neutrality is provided by monovalent 
negative ions). In this case one has to compare energies 
of two states. The first one consist of NaCl-like crystal 
consuming concentrations c_ of both positive and nega- 
tive spheres. The rest c_ positive spheres are free. Such 
state has the energy — 1.74i?c_ per unit volume of so- 
lution. The second state consists only of c_ triplets, in 
which two positive spheres are attached on opposite sides 
to a negative one. Energy of such a complex is —1.5E. 



We see that second state looses to the first. It can be 
shown that at any c+ /c_ , the largest possible fraction of 
spheres always self assembles in neutral NaCl-like crystal. 
Thus we are coming to conclusion that triplets (which 
are positive at c+ > c_ and negative at c+ < c_) do not 
form. Triplets can be considered as analogs of free and 
charged DNA-spheres necklaces. So in this case, neutral 
condensed state always dominates and there is no reen- 
trant condensation at all. 



VII. CONCLUSION 



In this paper, we presented a phenomenological theory 
of reentrant condensation. The theory is applicable not 
only to the solution of DNA and spheres but also to a 
much broader range of systems. 

Here, we would like to discuss possible implications 
of our theory for the natural chromatin assuming that 
the beads-on-a-string structure of 10 nm chromatin fiber 
is indeed determined by electrostatic interactions. In 
this case, the finite distance between nucleosomes tells us 
that the whole 10 nm fiber is charged. Electrophoretic 
experimentfS show that in low salt conditions, the net 
charge of 10 nm fiber is negative. (By the net charge we 
mean bare charges of DNA and histones plus charges 
which are Onsager-Manning condensed on them and, 
therefore, arc bound to the fiber with the energy larger 
than ksT). An additional argument for the negative sign 
of the charge is that the increasing salt concentration con- 
denses it into the 30 nm fiber. This means that 10 nm 
fiber obeys the scenario discussed in Sec. V for the point 
1 of Fig. |. 

Looking at the phase diagram of Fig. || we can ask 
whether one can increase concentrations of both octamers 
and salt much further so that that the point {s,p) be- 
comes higher than the decondensation curve Sd{p) and 
instead of 30 nm fiber or higher order structures we can 
get a positive stable 10 nm fiber. This is not easy be- 
cause free octamers decay into smaller histone dimers or 
tetramers. To keep equilibrium concentration of free oc- 
tamers s in solution at the necessary level, concentrations 
of histones should be large. In spite of this difficulty we 
suggest to try to create and study a positive 10 nm fiber. 
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